The conical point in the ferroelectric six- vertex model 



Dirk Jan Bukman and Joel D. Shore 

Department of Physics 
Simon Fraser University 
Burnaby, British Columbia 
Canada V5A 1S6 

February 6, 2008 



Abstract 

We examine the last unexplored regime of the asymmetric six-vertex model: 
the low-temperature phase of the so-called ferroelectric model. The original 
publication of the exact solution, by Sutherland, Yang, and Yang, and vari- 
ous derivations and reviews published afterwards, do not contain many details 
about this regime. We study the exact solution for this model, by numer- 
ical and analytical methods. In particular, we examine the behavior of the 
model in the vicinity of an unusual coexistence point that we call the "conical" 
point. This point corresponds to additional singularities in the free energy that 
were not discussed in the original solution. We show analytically that in this 
point many polarizations coexist, and that unusual scaling properties hold in 
its vicinity. 



Keywords: Six-vertex model, phase transitions, solid-on-solid models. 



1 Introduction 



The six-vertex model is one of the few exactly solved models in statistical mechan- 
ics. Its origin lies in the problem of counting the number of configurations of the 
hydrogen atoms in hydrogen-bonded crystals, such as ice ® and KH2P04.(@' i i In 
1967, Lieb, Yang, and Sutherland® obtained exact solutions for various versions of 
the model, and the most general solution was described by Sutherland, Yang, and 
Yang.® These publications contain virtually no details, and later more explicit ac- 
counts were published by Lieb and Wu,® and Nolden.*0' B The paper by Nolden 
contains a detailed description of the exact solution, and explicit results for some spe- 
cial cases where these can be obtained analytically. However, the low-temperature 
ferroelectric regime is not treated there. 

A more recent application of the six-vertex model is as a model of crystal surfaces. 
Various restricted solid-on-solid models can be mapped onto the six-vertex model, as 
was shown by van Beijeren for the body-ce„tereii;rtic lattice.® and by Jayapraka=h 
and Saam for the face-centered-cubic lattice. This last mapping in particular 
focuses attention on the low-temperature ferroelectric regime of the six-vertex model. 
In ref. ITUlit was conjectured that there would be a discontinuity of slope at a particular 
point in the crystal shape. Since in this mapping the crystal shape is, up to a 
scaling factor, given by the shape of the free energy surface as a function of applied 
fields, ^13) this implies the presence of a singularity in the free energy at this point. 
However, such a singularity was not mentioned in the general solution of the six- 
vertex model <B Our re-e..mn..tion of the exact solution m thi= regime, following 
the method presented in ref. shows that there is indeed a singularity at this point, 
and that it has new and unusual properties. In this point, a one-parameter family of 
polarizations coexists, and near it interesting scaling properties are found. 

The organization of the rest of the paper is as follows. In Section |^, we give a 
short, general discussion of the six- vertex model, and briefly review the expressions 
from the exact solution that we use further on. Then, in Section 0, we apply those 
expressions to a particular point in the phase diagram in the low-temperature regime, 
where a completely analytical solution is possible. The behavior of the free energy 
and the order parameters at and around this point is discussed. In Section ^ the 
same is done, in less detail, at the transition between the low- and high-temperature 
regimes. Conclusions and a summary of the results are given in Section ^. The 
appendices contain calculational details that were omitted from the main text of the 
paper. A brief account of this work was published in ref. |T^. 



2 The ferroelectric six-vertex model 



A configuration of the six- vertex model is given by a covering of the bonds of a square 
lattice with arrows, satisfying the ice rule: Every lattice point must have two incoming 
and two outgoing arrows. There are six possible vertices that satisfy this rule (see 
Fig. |l]). To make the model amenable to a solution by transfer matrix techniques 
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it is necessary to impose periodic boundary conditions in the horizontal direction. 
As a consequence, the number of down arrows is equal in each row of bonds; this 
conservation law is essential for the exact solution of the model. The transfer matrix 
method used in the exact solution is discussed in, e.g., refs. |, ^ |. 

After assigning an energy to each of the vertices, the partition function of the 
model is simply given by the sum over all allowed configurations of the Boltzmann 
factors of the configurations. The usual way to write the energies of the vertices is 
given in Fig. |l]. The quantities h and v are to be thought of as fields, acting on the 
horizontal and vertical arrows respectively, favoring one orientation over the other. 
(The name "asymmetric six-vertex model" refers to the fact that these fields create 
an asymmetry between left and right and between up and down arrows.) The fact 
that vertices 5 and 6 have the same energy is no restriction: since they are sinks 
and sources of horizontal arrows they have to occur in pairs because of the periodic 
boundary conditions, so only the sum of their energies is relevant. 

The variables that are conjugate to the fields h and v are the polarizations x 
and y. These are defined as a; = 1 — 2/^ and y = \ — 2/^, where, e.g., is the 
fraction of the horizontal arrows pointing to the left, x and y are order parameters 
for the various ordered phases. The free energy per vertex can be viewed as either a 
function of the fields, F{h,v), or of the polarizations, F{x,y). The two are related 
by a two-dimensional Legendre transformation 

F{h,v) = min{F{x,y) — hx — vy} (2-1) 

which implies F{h, v) = F{x, y) — hx — vy, with the fields given by 

^ ^ dF{x,y) ^ ^ dF{x,y) 
dx ^ dy 

or conversely 

y--'-^ (2.3) 

The free energy that is calculated in the exact solution is actually F{h,y), because 
of the conservation law for the number of down arrows. This amounts to solving the 
model for fixed values of h and y. The relevant free energies F{x,y) and F{h,v) are 
then obtained from F{h, y) through one- dimensional Legendre transformations. 

Most of the analysis in ref. |^ is valid for any values of the vertex energies. However, 
certain transformations that are applied do depend on them; more precisely, they 
depend on the parameter A = ^{rj + r]~^ — e^^'^), where rj = e^^ , and = l/fc^T is 
the inverse temperature. The discussion in ref. ^ is restricted to the case A < 1; we 
will focus on A > 1. There is a different transformation for A = 1 and for A > 1; 
both are given in ref. |^. As a consequence, various quantities and functions that occur 
in the exact solution have different forms in these two cases; they are given in table | 
(the expressions for the other ranges of A are given in ref. 0, 0). A > 1 corresponds to 
the low-temperature phase of models that are "ferroelectric", i.e. that have e < 5/2, 
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so that the vertices 1 and 2, with a non-zero net polarization, have lower energies 
(for h = V = 0) than the vertices 5 and 6 with zero net polarization (we choose 
5 > without loss of generality). By varying T at fixed S and e, a transition to a 
high-temperature phase (with A < 1) takes place at a critical temperature T = Tc 
(corresponding to A = 1). 

We now give a brief review of the main results in Nolden's paper® that we will 
use; for derivations and more details the reader should consult that paper. For fixed 
e, 6, and T, the exact solution provides expressions for the free energy, and various 
other quantities of interest, like the fields and the polarizations, as a function of two 
parameters, a and b. Different values of a and b then correspond to different points 
in the phase diagram, i.e. different values of h, v, x, and y. The first step is to find 
the function R{u) from the integral equation 

R{u) + ^ r K{u - v)R{v)dv = ^{u) (2.4) 

ZTC J -a 

R is a. function of u, a real parameter in the interval [—a, a]; it also depends on 
the values of a and This equation results from the consistency relations for the 
wavenumbers in the Bethe ansatz. Once R is known, the free energy per site can be 
obtained as the largest eigenvalue of the transfer matrix. 



- pF{h, y) = max 



±-(lnr/ + 2(3h) + — T (!>^^^{u)R{u)di 
2 2n J -a 



(2.5) 



with upper or lower sign for the right (R) or left (L) eigenvalue, respectively. Expres- 
sions for h and y can be obtained from the "generalized normalization" 

^(1 - y) - 2iph = p^{a) - 7^ r e(a - v)R{v)dv (2.6) 

2 271 J -a 

To find f , X, and the relevant free energies, F{h., y) must be differentiated with respect 
to h and y, on which variables it depends only implicitly through a and b. All this is 
described in ref. |^, and that treatment applies to the ferroelectric case without any 
modification. 

For general values of a and 6, the only way to solve these equations is numerically; 
the integral equation can be solved with standard techniques (See ref. ^ Chap. 18), 
and all quantities can be calculated. The ranges of a and b in table | lead to that part 
of the phase diagram where y > 0. The other half can be obtained through a simple 
symmetry operation. An analytic solution is possible for a = and a = — tt, and an 
expansion around these points can also be made. The case a = was examined in ref. 
TO; it corresponds to second order boundaries between completely and incompletely 



polarized phases in the phase diagram. We will focus on the case a = — vr , which 
corresponds to the two "conical" points. 



^Here and in what follows, explicit dependence of functions on a and b has been omitted for 
brevity. 
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3 Analytic solution for A > 1 



In the case that a = — vr, the integral equations in the previous section can be solved 
analytically by means of Fourier series, as was done for a = vr in the case A < — 1 in 
ref. 1^. It is then also possible to calculate an explicit expression for the free energy, 
and to examine its behavior around the point a = — tt. 

3.1 Finding R{u) 

The first step in the solution is to solve Eq. ( |2.4| ). To find the behavior both at and 
around a = —tt, we write the solution R as an expansion in[| e = vr + a: 



Riu) = Roiu) + J2 e^'SRM (3.1) 

m=l 

All integrals of the form /~" f{y)dv can be expanded in e as follows: 

f{v)dv = / f{v)dv 

+e J -IT 

oo 1 

+ E ^ {{-^r f'^-'\^) - e"^/("^-i)(-vr)} (3.2) 



m=l 



Using this in Eq. ( p.4| ), and substituting the expansion ( |3.1| ) everywhere, we get a 
set of integral equations, one for each order in e, containing Rq and the SRm- All 
functions Rq and SRm are defined on the interval u G [— tt,-/?], and all integrals are 
over the same interval. Since K and ^ have period 27r, the solutions Rq and 6Rm. 
also have this periodicity. (Note that by this periodicity terms with even m in the 
expansion ( |3.2| ) vanish. The terms with odd m reduce to twice the value at f = vr.) 
This periodicity makes it possible to solve for the functions Rq and SRm by applying 
a Fourier expansion to the equation. We will now examine the first five orders in e, 
and solve for _Ro and the first four 6Rm- 

The Fourier components of, e.g., Ro{u) are defined by 

oo 

Ro{u)= J2 ^ne-*"" (3-3) 

n=— oo 

and similarly for 6Rm, K and ^. The Fourier coefficients Kn are the same as those 
given in ref. ^ and ^„ is calculated in [Appendix A| : 



{n>0) 
2e'^ sinh(nz/) {n < 0) 
-2e^" sinh(nz/) (ra > 0) 
(n < 0) ^ ' 



^This e should not be confused with the energy e introduced earher. 
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Eq. to order e° is 

Roiu) r K{u - v)Roiv, b)dv = 

Fourier transforming this gives 

Rn[l-Kn]=in 

SO, for n 7^ 0, i?„ is given by 



(3.5) 
(3.6) 



Rr, 



° (b>u) 

_^{b-u)n ^^<o) ^^^""^ 
_^(b+u)n > 0) 

(ra < 0) 



{b < -v) 



(3.7) 



For n = 0, (p.6|) reduces to = 0, so -Rq is not yet determined. The Fourier series 
( p.3|) for i?o can be summed, and we find 



Rvsiu] = Ro- 



with the upper (lower) signs for b > u {b < —y)- Ro will be determined by the next 
order in e. 

Next, for we find 



SRAu) 



K{u - v)6Ri{v)dv = [K{u - 7r)i?o(7r)] 

TT 



271 J-n 

After Fourier transforming, and using K(u — n) = I]„(— l)"-ft'„e~*"", we have 

(-1)" 



(5i?i)„[l - K„ 



TT 



(3.9) 



(3.10) 



Since Ko = 1, there is a problem unless the r.h.s. is zero for n = 0. Thus Roi'n') must 
be zero. This requirement fixes Rq, giving 



Rn 



{b>u) 
{b < -u) 



(3.11) 



Once we put Roiii) = it is obvious that (5i?i)„ = for n 7^ 0. (5i?i)o will be 
determined in the next order. 

We continue in this spirit, finding expressions for the next terms: 



mi)o 

5R2{u) 
5Ri{u) 





6R2 



6 



ini-irK^R'oM /Q^ 
37r(l-i4-„) ^ 

{n = 0) 



(3.12) 
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The above expressions contain the first, second and fourth derivatives of Rq{v) 
evaluated at f = vr. These are obtained from Eq. (|3.^ ), 



±i 



1 + e^'^^Y 
(1 + e^Tfc)3 

(1 + e'^TMs 



(3.13) 



with the upper (lower) signs for b > u {b < —u). 
To recapitulate, the total expression for R{u) is 



R{u) = Ro{u) + e^5R2 + e'^5R^{u) + eHR^ + 0{e') (3.14) 
where 5R2 and 5R4, are real and independent of u, and SR^lu) is imaginary. 



3.2 Calculating y and h 

Having determined R, we can now determine y and h using (|2.6| ). We will first 
calculate the values yo and ho, the first terms in the expansion in e, and then find 
terms up to order from a slightly simplified equation. From table |, p^{—tt) is given 

by 



'i In 



+ 1 



The second term on the right hand side of ( |2.6| ) is 

L r e(-7r - v)Ro{v)dv = ^T.Rn r e(-7r - t;)e-*™dt; 
Zn J Zn „ J -TV 



(3.15) 



(3.16) 



The integrals /„ = — /^^ 0(— vr — v)e "^"-^dv are calculated in [Appendix B| . They turn 
out to be 



27r2 



:-l)"^(l-e 



-2nv\ 



2nu\ 



{n > 0) 
(n = 0) 
(n < 0) 



(3.17) 



For the aforementioned term we now have 

- 7^ E InRn = -vri?o ± ^[ln(l + e"^') - ln(l + e-"^')] 

with the upper (lower) signs for 6 > z/ (6 < —u). The final result then is 



(3.18) 



7r 



l-Vo)- 2iPho 



-ttRq + iv (b> u) 
-ttRq — w {b < —u) 



(3.19) 
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So we find, using Eq. ( p. Ill ), 



yo 

ho 



tanh(^) {b > v) 
-tanh(^) {h<-v) 

-v/213 {b > v) 
u/2[3 {b < -u) 



(3.20) 



Now we use Eq. (2.28) from ref. to find the next orders in e. This equation is 



- ^daV - 2i/3dah = R{a) - ^Q{2a)Ri-a) 



+ 



27T Ja 



6(a — v)daR{v)dv 



(3.21) 



Substituting the expansions for i?, daR, expanding all occurrences of a to order e'^, 
and using [Appendix B| , we find 



day - 2iPdah = e^R'^in) + e'^R'^irr) + .^e^^'V) 



TC 



30 



(3.22) 



Since Rq and R^^ are real and R'^ imaginary, we have 
day = 



leK{7r)-^e'Ri'\n) + 0{e') 



2i3dah = ^e'?^ + 0{e'') 



(3.23) 



Putting all this together, we now have power series in e for y and h: 



with 



y 

h 



yo 
y2 

2/4 

ho 



yo + e%2 + e%, + Oie') 
ho + e'hs + 0{e') 



(3.24) 



± tanh 



'b =F i^^ 



1 e^^^l - e^^^) 
3 {l + e'^^^y 

1 e^^^(l - lle^^^ + lle2('^Tf>) _ e3(^T6)) 
~60 

V 

1 P^^^ 



(1 + e^^^^f 



67r/3 (1 + e^T^)2 
with the upper (lower) signs ioi b > v {b < —u) 



(3.25) 
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3.3 Calculating x and v 

To rigorously calculate expressions for x and v, and the Legendre transforms F{x, y) 
and F{h,v), one has to construct the derivatives of F{h,y) by applying the chain 
rule. This is what is done in ref. |^ However, by making a few assumptions, it is 
possible to derive the same expressions from the symmetry properties of the exact 
solution. The assumptions can be checked against known exact results, in which case 
they turn out to hold. 

The first step is to realize that, instead of calculating the free energy as a function 
of h and y, one might just as well calculate it as a function of another pair of variables, 
for example v and x. The results will be exactly the same as before, except for the 
change of variables, and the region of the phase diagram that one finds the solution 
for. The preceding calculation gives h, y, and F{h, y) expressed in terms of the 
parameters a and b. The solution holds for y > 0, which corresponds to regions II 
and III (for b > u), and region IV (for b < —v) in Fig. 0.0 The solution for y is of 
the form 

_ j Yi(a, 6) in regions II and III , . 

\ ^2(0, &) in region IV ^ ' ^ 

where Yi and Y2 are given functions. If, on the other hand, one calculates F(f,x), 
one finds exactly the same solution, but now v and x play the roles of h and y, and 
the parameters are now called a and b to distinguish them from the other set of 
parameters. The solution is now for that part of the phase diagram where a; > 0, 
regions V and IV (for b > u), and region III (for b < —v). The expression for x is 
just the same as the one found previously for 

_ J Fi(a, b) in regions V and IV , . 

[ Y2{a,b) in region III 

with the same functions Yi and Y2. We now have expressions for both y and x in 
regions III and IV, only they are given in terms of two different sets of parameters, 

a, b and a, b. We need a relation between them to find x and y expressed in a and 

b. A simple ansatz is to assume that a = a, and that b and b are related linearly. 
Since the line b = (pQ corresponds to 6 = —u, and b = 00 to b = —00, this would give 
b = —b + (po — u. Substituting this, we find for x 

^ ^ i Xi(a, b) = Y2{a, -b + (po - v) in region III ,^ 

I X2(a, b) = Yi{a, —b + (po — u) in region IV \ ■ ) 

If we further assume that the expression for x in region III also holds in region II, 
as is true for y, we now have expressions for x throughout the region y > 0. The 

^In this section we are only concerned with the so-called incompletely polarized part of the phase 
diagram, which consists of the regions I - VI. This is because this is the only part that is mapped 
out by varying a and h. In the rest of the phase diagram, x and y are frozen at their extreme values 



±1. This is discussed in detail in section 3.5 
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expression for x for h > u can simply be obtained by taking y for h < —v and 
substituting —h^(\)Q — v for h. Similarly, x for b < —u follows from the expression for 
y for h> V with the same substitution. The same prescription can be used to obtain 
V from h. 

It is possible, to a certain extent, to check the assumptions we made. As a first 
check, we have calculated x and v to zeroth order in e the rigorous way, obtaining 
the same results as above. Second, from the symmetry properties we can also derive 
consistency relations for the expressions for y in various regions. For example, we can 
repeat the above reasoning for F{—v, —x). This gives the region of the phase diagram 
where x < 0, regions I and II (for h > v), and VI (for h < —v). The expressions for 
X are now 

^ ^ f -Fi(a,6) in regions I and II ,^ 
1 —Y2{ci^h) in region VI 

If we make the same assumptions as before, a = a and h = u + (po — b. Then in 
region II, x is given by —Yi{a, —h + v + 0o)- But in that region x is also given 
by 1^2(0, + 00 — '^)- So to be consistent the expressions for y have to satisfy 
Yi{a, —h + v + 0o) = — ^2(ct) —h + (f)Q~v). Checking the expressions for the first 
few orders in e ( p. 251 ) shows that this is indeed true, giving confidence that the 
assumptions we made are justified. In addition, the numerical results confirm the 
expressions for x and v found in this way. 



3.4 Finding the free energy 



Now we are in a position to calculate F{h,y) from ( |2.5| ). We need to know the 
integrals 



du 



These can be expressed in terms of integrals of the form 



J n 



dz 



O 



In 



z-A 
z-B 



+ C 



(3.30) 



(3.31) 



for n > 0. These are calculated in [Appendix (J 



The details of the calculation of the zeroth order part Fq of the free energy as a 



function of y and h can be found in [Appendix D| . There it is shown that the right 
eigenvalue is largest for 6 < —u and b > (po, and the left eigenvalue is largest for 
u < h < (pQ. Even though at 6 = 0o the two eigenvalues cross, the analytic form of 
the expression for the free energy Fq is the same on either side of 6 = 0o- (In fact. 



this remains true at least up to order e ). In [Appendix D[ we find 



(3.32) 
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with the upper (lower) signs for h> v [h < —v). Note that for b > u, vq = dF/dy = 
so a = —71 and b > u corresponds to the point {ho,Vo) = (— Z//2/9, Z//2/3). 
Similarly, a = — vr and b < —u corresponds to (/i^'^o) = i.^l'^f^) ~^^/2/3). 

To find the next two orders of F{h, y) in e we expand Eq. ( p.5| ) as usual, to order 
e^, which gives 

- - F„) = ± ii ; ; fz 

with the upper (lower) signs for 6 > z/ (6 < — z/), and a = b — (pQ + u. The general 
structure of F{h, y) is thus found to be 

F{h, y) = Fo + e^F2 + e'F^ + 0{e') (3.34) 

The coefficients F„ can be expressed in terms of the coefficients in the expansions of 
h,v,x, and y (see Eq. ( p.25| )): 

6 

Fq = -- + voVo 
F2 = voy2 

F3 = -hxo (3.35) 

Since the Legendre transforms of the free energy are given by F{x, y) = F{h, y) + xh 
and F{h, v) = F{h, y) — vy, we also have 

F{x, y) = -- + voyo + hoXo + (fo2/2 + hoX2)e^ + O(e^) (3.36) 
Fih,v) = -^--{hxo + v,yo)e' + 0{e^) (3.37) 

3.5 Discussion 

The overall phase diagram as a function of h and v can be obtained numerically; its 
central region is shown as the base plane in Fig. ^. It contains four phases where 
the polarizations {x,y) are "frozen in" at their extremal values (±1,±1). Two of 
these phases, (1, 1) and (—1, —1), meet along the line h = —v. x and y are discon- 
tinuous across this line, so it represents a first order phase transition. In between 
the four frozen regions are two incompletely polarized phases, where x and y change 
continuously with h and v. They are separated from the frozen phases by continuous 
transitions, which have been shown to be of the Pokrovsky-Talapov (PT) type. 
The points where two PT lines meet the first order line are the "conical" points. 

To show the relation between the values of the parameters a and b, and the 
phase diagram in terms of h and v, several lines of constant a and b have been 
drawn in Fig. ^ This relation is shown only for the part of the phase diagram with 
b > {v + (po) /2. This gives one quarter of the phase diagram; the rest of the diagram 
is then obtained by symmetry. We restricted our computation to this region to avoid 
certain numerical difficulties that arise near b = ±z/: At 6 = ±i/, R{u) has a pole and 
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this makes the solution of the integral equations increasingly difficult a.s b ±z/. An 
additional difficulty is that at 6 = (pQ, $'^'^(-u) have poles. This does not affect the 
integral equations themselves, but does make the integrals involved in the calculation 
of F{h, y), X, and v increasingly difficult as 6 — *• 0o. We have dealt with this problem 
by using an integration routine with adaptive stepsize control (See ref. |l^, Chap. 4), 
as opposed to the simple Gaussian quadrature used for solving the integral equations 
and determining the less-sensitive integrals. However, near b = (po the integration 
does become quite slow. 

Fig. ^ also shows the free energy surface F{h,v) superimposed on the phase di- 
agram. F{h,v) is purely linear in h and v in the frozen phases, reflecting the fact 
that {x,y) = —{dF/dh,dF/dv) is fixed. In the incompletely polarized phases, x and 
y change continuously, so F{h,v) is smoothly curved. Along the PT lines the curved 
part of F joins smoothly onto the linear part, since the singularity in F{h, v) corre- 
sponding to the (second order) PT transition leaves the ffist derivative continuous. 
Along the ffist order transition line there is a jump in slope, so that a ridge is formed. 
At the two conical points the slope of F depends on the way in which the point is 
approached, as will be shown later on. 

We now turn to the phase diagram in x-y space. Fig. |^. Here the four frozen phases 
correspond to the corner points (x, y) = (±1, ±1). The incompletely polarized phases 
take up the rest of the phase diagram outside the lines and The two conical 
points now correspond to the two lines and which are given by (xo(&), ?/o(^)) 
for b > V and b < —u, respectively (see Eq. ( |3.25|) ). By eliminating b in Eq. ( p.25|) 
we find that the lines i± are given by 

xo ±tanh[(0o - z^)/2] 

^0 = TT r~u72 ^^-^^^ 

1 ± Xo tanh[(0o — '^j/^J 

From Eq. ( p.36| ) we see that on, e.g., where a = — tt and b > u, the free energy 
is F{x,y) = —6/2 + v^y^ + HqXq = —6/2 + (z//2/3)(?/o — Xq). This means that the 
slope of F{x,y), which is equal to {h,v) = (— z//2/?, Z//2/3), is the same everywhere 
along £+, and that every (xo,i/o) on corresponds to the same point in {h,v), the 
conical point. So, all the polarizations (xo,yo) are stable in the conical point. This 
is similar to the coexistence of a one-parameter family of magnetizations m with 
|m| = mo{T) in a zero-field XY model for T < Tg. However, there the coexistence 
is due to an obvious symmetry in the Hamiltonian, while in the six-vertex model the 
symmetry that allows the polarizations to coexist is not an obvious symmetry in the 
Hamiltonian and therefore appears to be generated spontaneously. 

It can be shown that different values of b correspond to entering the conical point 
with different angles from the incompletely polarized region; if b is kept fixed as 
a — > — TT, the conical point is approached from an angle 

cosh(6 — z/) — cosh(6 — 0o) 



9 = arctan 



(3.39) 



2 -|- cosh(6 — u) + cosh(fe — 0o 

with respect to the line h = —v. Since every different b leads to a different value of 
(x, y), and (x, y) gives the slope of F{h, v), the slope in the conical point depends on 
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the angle at which it is entered. This gives F{h, v) a geometry similar to the tip of 
a cone at this point, hence the name "conical" point. 

The exact solution does not give any states with polarizations in the regions 
between and Thus, no pure equihbrium states with those polarizations exist, 
and strictly speaking the free energy F{x^ y) is not defined in this region. However, 
a state with an average polarization (a;, y) in, e.g., can be formed as a mixture of 
states on £+, properly weighted to give the right average polarization. This leads to 
a free energy that is also an average of the free energies of the states in the mixture, 
which means that it is linear in x and y. In this way, the regions are interpreted 
as describing mixtures of the coexisting phases on similar to states that are a 
mixture of coexisting gas and hquid phases at a gas-liquid phase transition. The 
slope of the free energy surface F{x,y) is then given by {h,v) = (— Z//2/3, z//2/?) in 
C+, and {h,v) = (z//2/3, — 1//2/3) in C~. This matches the slopes on the lines 
and so that the first derivative of the free energy F{x,y) is continuous across 
these lines. Its higher derivatives are not continuous, however, and the free energy 
has singularities at the lines i±, since an analytic function cannot be matched to a 
purely linear function in a regular way. This form of the free energy also means that 
the fields corresponding to the regions are the same as those for the lines i±, i.e. 
the mixtures described by occur in the conical points. The two linear parts of 
F{x, y) meet along the line x = y, forming a groove. This groove corresponds to the 
coexistence of the two frozen phases along the first order transition in Fig. |^, and it 
was already mentioned in ref. ^. However, it is stated there that the groove is the 
only set of points in x,y & (—1, 1) where F{x,y) is non-analytic, while in fact there 
are additional non-analyticities along the lines i±. 

To examine the nature of the singularities in the free energy we now specialize to 
the case 6 = (z/ + 0o)/2, so that we are on the lines y = ~x and v = —h, moving from 
the incompletely polarized phase towards the conical point as e ^ 0. The behavior 
along this line should be qualitatively the same as along any other line emanating 
from the conical point; however, it is much easier to deal with just two variables 
instead of four. Solving for e in Sv = v — vq = e^v^ + O(e^), and substituting this in 
the expression for y (Eq. p.25|) , we find that 



yi^) = yo + ^Sv'/' + 0{6v'/') (3.40) 



V 



3 



So, as 6v goes to zero, y approaches its value of yo in the conical point with a slope 
that diverges like 6v~^^^. This value is in agreement with recent work*^E' H' which 
shows that certain (l+l)-dimensional surface growth models map onto the conical 
point in the {h,v) phase diagram, and that Kardar-Parisi-Zhang (KPZ)(0^ scaling 
holds. In particular, the exponent of 2/3 in Eq. ( |3.4CI| ) is given by where 
z = 3/2 is the dynamic exponent for the KPZ universality class. Relation (|3.40| ) can 



be integrated to give F{v), since dF{v)/dv = —2y (the factor 2 comes from the fact 
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that X = —y and h = — f). This gives 

F{v) = - 2yo6v - -^6v'/' + 0{6v'/') (3.41) 

The first few terms of this expression reproduce the expansion of F{v) given in terms 
of e earher (Eq. |3.37|) . Similarly, an expression for F{y) can be derived by writing v 
in terms oi 5y = y — y^ and integrating the relation dF[y)/dy = 2v. It is 

Fiy) = -^- + 2vm + 'ivoSy + -^,5y''^ + 0{5y"^) (3.42) 



Again, this confirms the expansion given earlier (Eq. |3.36 ). 



The appearance of the terms 5v^^^ and Sy^^"^ means that both F{x^ y) and F{h^ v) 
have non-analyticities at the first-order transition, in addition to the expected discon- 
tinuities in their derivatives. Thus, the traditional mean-field picture of a first-order 
phase transition as two branches of the free energy simply crossing is not valid. As 
the transition is approached from the incompletely polarized phase, the free energy 
becomes singular, and there is no obvious way to extend this branch of the free en- 
ergy into the regions in Fig. ^ to find metastable states. It is thought that the 
existence of non-analyticities at coexistence boundaries is a generic feature; the 
two-dimensional Ising model, e.g., shows a weak essential singularity in F{h) at /i = 
and T < T^. Vector spin models (where, as here, a whole continuous family of order 
parameter values coexists at a single point) show power law singularities similar to 
the ones found here, albeit with different values for the exponents (we are not aware 
of any previous observation of, e.g., a "susceptibility" d'^F{v)/dv'^ with an exponent 
of —1/3). As a further comparison we mention that at the PT transition, where y 
goes to 1, the free energies behave like^© 

Fiy) = A + B{l-y) + Cil-yf + --- 

F{v) = a + P6v^^^ + -f6v^ + ■ ■ ■ (3.43) 

Note that the conical points are reached either by letting a — vr or by letting 
b —>■ ±oo (or some combination of both), as was remarked earlier. Taking a —>■ — tt 
brings us into the point along some nonzero angle with respect to the PT line in {h, v) 
space, whereas b ±oo takes us in tangent to the PT line. The scaling properties 
derived above [Eqs. ( |3.40| - P^^2D ] correspond to the former case. In the special case 



of a trajectory tangent to the PT line, the scaling would be determined by the limit 
b ±cxo. In particular, we know that the free energy F{h,v) is analytic in h and 
V along the PT line itself (6 ^ oo with a = 0). By coming into the conical point 
along other curves which are tangent to the PT line, intermediate scaling forms of 
various kinds could be obtained, but such scaling behavior is not generic and therefore 
somewhat contrived. 

This brings us to a closely-related point which is: In which region is the behav- 
ior near the conical point dominated by the proximity to the conical point, and in 
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which region is it dominated by the proximity to the PT hnes? For example, we 
might ask how to determine the region over which we can observe the scahng {e.g. 
the Pokrovsky-Talapov behavior) associated with the boundary of the completely 
polarized region. The preceding discussion suggests that this region is very narrow. 
Indeed, by considering when the second term in the denominator of (A19) of ref. |10| is 
of order 1, we find that a curve of the form df ~ separates the two regions. Here, 
dc is the distance in the h-v plane to the conical point and df is the perpendicular 
distance to the PT line. Note that this curve is tangent to the PT line, deviating 
only quadratically as one moves away from the conical point. Application of (A19) 
in ref. |1^ also tells us that the coefficients of the {1 — y) and (1 — y)^ terms of the 
expansion of the free energy about the completely polarized region [as in Eq. ( |3.43| )] 
are unchanged by the presence of the conical point, and thus that the coefficient of 
the PT scaling term is not modified. However, higher-order terms in the expansion 
diverge as one approaches the conical point, representing the fact the the PT scaling 
will be observable over an ever-smaller angular region near the PT line. 

We can also examine the temperature dependence of various quantities. For ex- 
ample, the jump in y along h = —v in the conical point. Ay = y^ip = (0o + '^)/2) = 
tanh((0o — ^)/^) behaves as follows: For T ^ 0, ^.e. A ^ oo, we find that 0o — ~ 
p{5 - 2e) oo, so that Ay ^ 1 - 2e-'^(^-2£)/2_ j. _^ ^ ^ 

00 - oc z/, and z/ oc A^/^ ^ (j^ _ ^-^at Ay oc (T^ - Ty/^. Qualitatively 

the same temperature dependences are found for the opening angle 26q between the 
two boundaries of the incompletely polarized phase at the conical point; can be 
obtained from Eq. (|3.39|) by letting 6 — oo, which gives tan 6^0 = tanh((0o — ^)/'^)- 

Finally, we might ask whether there are any analytic approximations we can make 
which are valid over the entire phase diagram. Indeed, we ca n pe rform an expansion of 
Eqs. ( |2.4| - |2.6| ) to lowest nontrivial order in the temperature, fi^-* To do this, we ignore 
terms in .^(m), K{u — v), p^{u), and which are down by factors of order e~'^'^ 

(or higher powers) from the dominant terms. Since to this order the kernel K{u — v) 
is identically 1, the integral equation can then be solved analytically and we can write 
down explicit expressions for h, v, x, and y in terms of a and b. The expressions for h 
and V, however, involve an integral which apparently must be evaluated numerically. 
Furthermore, even the expressions for x and y cannot be inverted to eliminate the 
parameters a and b from the equations entirely (leaving us, as we had hoped, with 
direct relations between h, v, x, and y). Therefore, these expressions, while providing 
us with a simpler numerical task than solving integral equations, do not lead to any 
profound simplification. 

Since we have ignored terms down by orders of e'"^", the low-temperature ap- 
proximation is quite accurate even up to a large fraction of the critical temperature. 
One particularly interesting fact about the accuracy of this approximation should be 
noted: Although the general features, such as the locations of the PT lines (corre- 
sponding to a — s> 0) are not obtained exactly, everything involving the conical points 
{i.e., the limit a — vr and the limit b ±oo) is exactly correct. This includes 
the location of these points in the {h,v) plane, the expressions for l±, and even the 



14 



expansions for x, /i, and v about the conical point [Eqs. ( |3.24| - py^ )] to the orders 
we know them. (In particular, this means that the transition temperature is ob- 
tained exactly — a curious fact for the first term of a low-temperature expansion!) We 
hypothesize that the fact that the behavior at the conical points is obtained exactly 
could be evidence that this behavior is controlled by a zero-temperature fixed point. 



4 Analytic results for A = 1 



The case A = 1 corresponds to the critical point T = Tc, where the first-order line 
in Fig. 1^ has shrunk to zero length (see Fig. and the two conical points at the 
end of this line have now coalesced into one point, which will disappear for A < 1. 
Here, it is possible to obtain some analytic results for a —>■ — oo, which corresponds 
to this limiting point. It turns out to be much harder to make an expansion about 
a = — oo than about a = —it, and we were only able to obtain a limited amount of 
information about this case. 

We can go through the same steps as in section ^ using the relevant expressions 
from table |I[ Now equation (|2.4| ) can be solved by means of a Fourier transform, 
defined as 



Ro{t) 



271 J- 



e''''Roiu)du 



(4.1) 



The Fourier transform of K is again the same as in ref. |^, and the one for ^ can easily 
be calculated: 



m 
m 



\t\ 



(t > 0) 

2e^* sinh(t/2) [t < 0) 
-2e^*sinh(t/2) (t > 0) 
(i < 0) 



(6 > 1/2) 
(b < -1/2) 



(4.2) 



Solving 



Ro{t) 



for -Ro(^) gives 








(t>0) 



.g(b-l/2)t (^<o) 



(t<0) 



(b > 1/2) 
{b<-l/2) 



(4.3) 



As usual, -Ro(O) is not determined. However, that does not prevent us from Fourier 
transforming back to find Ro{u), 



i?n(u) 



(■oo 



b— 2 —iu 



b+^—iu 



ib<-i: 



(4.4) 
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It is, to a certain extent, possible to expand R{u) in 1/a around a = — oo. It can 
be shown that the leading corrections are of the form 

R(u) = Ro(u) H h I \ (4.5) 

a 

with / and g real functions, / even and g odd. These two functions are determined by 
two integral equations which still contain a. Even though numerical evidence shows 
that / and g have well-defined limits for a — — oo, there is no obvious way to take 
this limit in the integral equations, and we have not been able to obtain closed-form 
expressions for the two functions. 

To get Uq and we again apply ( |2.6| ). It turns out that we have to be slightly 
careful in taking the limit a — > — oo, because we have to consider the possibility that 
h diverges along with a, so we cannot just set a = — oo as we did before. From table | 
it is apparent that p^{a) — > no matter how a and b go to infinity. For the integral 
term we find 

. „ , . . 1 



— / Q{a - v)Ro{v)dv ^ — 0(-oo) / Ro{v)dv 



= arctan(^^^) (4.6) 

with the upper (lower) signs for b > 1/2 (6 < —1/2). So we find that for a — oo 

_\ l-l arctan(-a/26) {b > 1/2) , . 

~ \ 1 - I arctan(a/26) {b < -1/2) ^ ^ 

and /lo = in both cases. So by letting b go to infinity along with a, any positive 
value of yo can be selected. It follows from the expansion ( [4.5|) that the corrections 
to both h and y will be of order 1/a; however, the coefficients of these terms cannot 
be calculated without knowledge of the functions / and g. 

To find expressions for x and v we can again use the symmetries discussed in 



section 3^, if every occurrence of is replaced by 1/2. In particular, we immediately 
find that Vq = 0. 

It is also possible to find the free energy in the limit that a — » — oo. As before, 
it turns out that for b < —1/2 and b > (po the right eigenvalue is largest, and for 
1/2 < 6 < 00 the left eigenvalue is largest. In all cases we find that Fo{h,y) = —S/2, 
independent of b. This confirms the above conclusion that vq = dFo/dyo = 0. 

We can again examine the scaling behavior of the free energies. As in the case 
A > 1, F{h,y) is a regular power series in the expansion variable, 1/a. However, in 
contrast to that case, the expansions for y and h both contain linear terms in the 
expansion variable. This implies that the polarizations are non-singular functions of 
the fields. Doing the Legendre transform of F{h,y) to obtain F{x,y), we find that 
the coexistence regions have shrunk to the line x = y, and the behavior around 
this line is now quadratic in x — y. So both the groove at x = y and the singularities 
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at the lines C.± (which are now collapsed to x = y) have disappeared. The free energy 
can thus be written 

y) = ~\ + A{x - yf + 0{{x - yf) (4.8) 

where the coefficient A = A[x + y) is a function of the position along x = y. Similarly, 
the free energy as a function of the fields is given by 

Fih, v) = -^- Voih + v) + Oiiv - hf, (v + hf) (4.9) 

where yo depends on the ratio a/b according to Eq. ( [4.7|) , i.e., it depends on the 
details of how the point h = v = is approached. Note that both the jump in the 
slope of F[h,v) on entering this point, and the singularities on approaching it, have 
disappeared; however, the slope in the h + v direction (given by —yo) still depends on 
the way in which this point is approached. This reflects the fact that, as T — > from 
above, the two PT boundaries approach each other, and the curvature of F{h, v) in 
the h + v direction goes to infinity. 

5 Conclusion 

To summarize, we have studied the phase diagram of the asymmetric, ferroelectric 
six-vertex model in the low-temperature phase. From numerical results we get a 
global picture of the behavior of the free energy as a function of the polarizations 
and fields, F{x,y) and F{h,v), and by analytical methods we found several new and 
interesting features. For T < Tc (A > 1), in addition to the groove that was already 
known to exist in the free energy surface F{x,y), there are additional singularities 
along the lines £±, corresponding to two conical points in F{h,v). In these points, 
which lie at the ends of a ridge in F{h, v), all polarizations on i± are stable for one 
value of the fields h and v. Also associated with these points are the two coexistence 
regions C^, where mixtures of the coexisting states are stable. 

As T —> Tc (A ^ 1), the size of the regions decreases, as does the length of 
the ridge in F{h, v). At T = Tc, the two conical points merge as the ridge is reduced 
to a point, and the regions vanish. There is still a family of polarizations, x = y, 
that is stable in the (one) conical point. If there are any non-analyticities left in 
the free energy they are of higher order than the terms we have been able to study. 
Above Tc, the conical point disappears completely. 

The behavior of the free energies when approaching the conical points is rather 
unexpected. Both the fact that there are power law divergences, while the phase 
transition is discontinuous and thus ffist-order, and the specific exponents of the 
power laws, are interesting. It would be useful to find metastable states in the regions 

with a uniform polarization, and examine how such a state decays into a mixture 
of stable states on the lines i±. However, because of the singularities in F there is 
no obvious way to analytically continue the free energy into the unstable region, in 
order to find such metastable states in a straightforward manner. Also, there does 
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not seem to be any mechanism to select one of the many possible mixtures that such 
a state might break up into. Thus one can ask how the process of phase separation, 
or spinodal decomposition, would actually take place. 

A particular system in which this issue needs to be addressed is the phase sepa- 
ration of crystal surfaces. ^^3' Through the mapping of the six- vertex model onto 
restricted solid-on-solid mo delsi' 0) the polarizations x and y are found to be equiv- 
alent to surface orientations, and the free energy F{h, v) turns out to be a replica of 
the crystal shape itself, C:^' c3 while F{x,y) is the surface tension. Thus, the conical 
points represent actual points on the crystal surface that have a conical geometry.^ 
In these points, various surface orientations all coexist. Also, there is a discontinuity 
in orientation as one crosses a conical point, since the orientations that correspond to 
states in the regions are not stable anywhere. Furthermore, if a crystal is cleaved 
to form a surface having one of these unstable orientations (say in C^), it should 
decompose into a mixture of stable orientations selected from the set £+, which then 
coarsens over time. This is the phase separation process mentioned in the previous 



paragraph. These issues are discussed in more detail in ref. 21 
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Appendix A The Fourier coefficients of ^{u) 
The Fourier coefficients of ^ are given by 

= ^ /- MMe"" = i r in (A.l) 

27r J-TT In J-Tx cosh u — cos(m + lo) 

For n > 0, substituting z = e^^ gives 

- .sinh ^ /" ^ ^" (A 2) 

71 Jq 2z cosh u — — z'^e~'^ 

The contour is a counterclockwise circle with radius 1. The poles are at zi = e~'^^^ 
and Z2 = e^^^ . For b > u both poles lie outside the contour, and ^„ = 0. For b < —u 
they are both inside the contour, and the integral is 



27ri- 



-i sinh V \ z^ 



TT 



+ 



2 sinh —2 sinh 



26"" sinh ni/ (A.3) 



''For A < — 1, the free energy as a function of x and y has a conical point {i.e. there is a cusp 
in the surface tension); this leads to a flat region in the crystal shape F{h,v): a facet. '■0 & This 
situation is thus the converse of what takes place for A > 1. 
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So for n> 



{b> u) 

— 26^"^ sinh(nz/) (b < —u) 



(A.4) 



For < 0, we substitute z = e ™ instead, and in a similar way we find 
f 2e''" sinh(nz/) (6 > u) 







(6 < -u) 



(A.5) 



Appendix B The integrals 

The integrals /„ are defined as 

J — IT 

The function Qi—Ti — v) can be rewritten as 



0(— TT — t)) = « < In 



+ 2z/ 



(B.l) 



(B.2) 



where the In has a branch cut along [0, oo), so that its imaginary part runs from to 
27rz. This is because 0(— tt — v) must run from for v = — vr to — 27r for v = tt. 
For n < we use ( |B.2|) , define z = — e™, so that 



-i-ir ! dzhn- — ^ + 
Jo I B - z 



(B.3) 



with A 



'2u 



B = e , C = 2v. The contour is given in Fig. M. So we must examine 



integrals of the type 
I 



n- , An^-^ + C[z"-'dz 
O I B - z 



[B.4) 



for n > 0. For n ^ 



In + 



+ 



X - A 



B-x 

, X- A 
In 



+ C^ + 7i 



B-x 



+ C + 72 



x'^-^dx + little loop = (B.5) 



where now the In's are the principal branch {i.e. with a branch cut along (— cxd,0]), 
the 7's are the imaginary terms that result from approaching the cut from above or 
below. Since the integral around the little loop goes to zero, and the line integrals 
give (1 — A") (71 — 72)/"",, we find 



27ri 



n 



[l-A-) 



(B.6) 
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For n = there is an extra pole at z = 0, so we must replace the right hand side of 



( pip with 



2ni{\n — + C) = 27ri(ln A-\nB + in + C) (B.7) 

The two line integrals give —(71 — 72) lnA = 27rz In A. So now 

Jo = -2n^ + 2ni{C - In B) (B.8) 

For < a similar calculation is performed by defining z = — e~™. The end 
result for J„ is 

- (n > 0) 

27r2 " ln = 0) (B.9) 

(-l)'^^(l-e2"'') (r2<0) 

Appendix C The integrals Jn 



Wanted are integrals of the form 



Jn= dz 

Jo 



1 ^ ^ ^ 
In - + C 

z — B 



z''-^ (C.l) 



where n > 0, A, i? > 0, and the In is the principal branch. There are six different 
cases, corresponding to the six orderings of A, B, and 1. They can all be computed in 
a similar spirit as in [Appendix B| . As an example, the contour for the case B < A < 1 



is shown in Fig. As another example, the contour from [Appendix B| can be used 
for the case A < 1 < B, provided that the values and 27r for the imaginary part 
of the logarithm along the branch cut are replaced by — tt and tt, respectively. The 
results for J„ are displayed in table fl[ 

Appendix D The zeroth order free energy 

In this appendix we calculate the free energy ( p.5[) , to zeroth order in e, for both the 
left and right eigenvalue, and compare them to find the larger eigenvalue {i.e. lower 
free energy). There are four cases: 

Case 1: 6 > z/, $ = $^ 

In this case, Rq is given by (see Eqs. ( p.3[) and ( p.7[) ) 

-1 

Roiu) =Ro- e(''-")"e-'™ (D.l) 
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Substituting this in Eq. ( ^.51 ) gives 

1 



E j_^^''iu)e-'''''duRn (D.2) 



Now, defining 2; = e*", and using [Appendix C 



^^(u)e-'"''du = -i - v \ z-^'-^dz 

^ ^ Jo[z- e-'?^o+^ J 

{b > 0o) 

_^j^g_„(;,_0o) _ ^_2^ + 00 < 6 < 0o) < 0) 

27rz/ (6 > 0o) 

27r(z/ - 00 + (-2z/ + 00 < & < 0o) = 0) (D.3) 

-27rz/ (b < -2z/ + 0o) 

Substituting ( [TISD into ([Hg) gives 

a) z/ < 6 < — 2z/ + 00 

-PFR{h,y) = ^{Inri +2ph) + uRq -\n{l - e""^") +\n{l - e^"-'^") (D.4) 

b) -2u + 00 < & < 00 

- = i(lnr/ + 2(3h) + (0o - - 6)i?o - ln(l - e'^"*") + 

ln(l-e'^-^) (D.5) 

c) & > 00 

- PFR{h,y) = ^{\nr] + 2[3h) - pR^ (D.6) 

Case 2: 6 > z/, $ = $^ 

A very similar calculation gives 

a) z/ < 6 < 00 



--(Inr/ + 2/3/i) - z/^o - ln(l - e^-"^") + ln(l - e-""^'') (D.7) 
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h) (f)o < b < 2iy + (po 



-PFL{h,y) = -^(lnr; + 2/3/i)-(0o + ^/-6)i?o-ln(l-e'^-^)+ln(l-e-'^-'^°)(D.8) 
c) & > 2i/ + 00 

- PFL{h,y) = ~{\n7] + 2(5h) + uRo (D.9) 

Case 3: 6 < -z/, $ = $^ 
Here we find that 

- (5FR{h, y) = ^(ln77 + 2(5h) + vR^ (D.IO) 

Case 4: 6 < -i/, $ = 

Here, 

- /3Fi(/i, y) = -^(ln77 + 2/3/i) - i/^o (D.ll) 



Now we have to collect all information on — /3FR^i(/i, y) and find the maximum over 
i?, L for the various intervals. The simplest case is 6 < —v: Here, 

- l3FR{h, y) = ^(Inr; + 2/?/i) + vR^ = I3Fl{K v) (D.12) 

Also, 

^o = -t4^ 2/?/i = i/ (D.13) 

X I 6 

We also know that luTy > i/, say \nr] — v -\- C, with C > 0. Then 

-/?F^(/.,y) = ^^-^ + ic>0 (D.14) 
So —PFfi{h, y) corresponds to the larger eigenvalue in this regime. Then, 

- m{K y) = Y + K^o + ^) = Y + (D-15) 

For h> v similar reasoning shows that for h G [i^, 4>q] the left eigenvalue is largest, 
and for b > (po the right eigenvalue is largest, and we find 

- pFoih, y) = f - K^o + ^) = Y " 
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Tables 



A = 

gip"(n) . 

Range 
Range 
0(m- 
Range 
K{u- 

00 : 



of a 
of 6 

v) = 

of e 

v) = 



A = 1 



-l+2m-26 

[-00,0] ^ 

'2' 



-OO, 



u 



2 arctan(u — v) 
(-7r,7r) 

2 

4 



l+4(u+fi6)2 
1 



^0 



In 



2 r;-l 
Q — 1+m— fc 

(/)0+m— fc 



A > 1 

cosh u (z/ > 0) 



[-7r,0] 

(— OO, — z/) U (z/, oo) 
2 arctan 



-27r,27r) 

sinh 2v 



cothz/tan 



cosh 2i/— cos(ii— d) 

sinh u 

cosh 1^— cos(n+i6) 

rj—e'^ 

-iu_^ — 2v + 4)Q—b 

Q — iu q4'0~^ 

-iu ^2u-\-4>Q — b 



+ z/ 
— u 



Table I. The definitions of various quantities and functions for the two transformations 
for A > 1. Note that some of the functions depend on b, but the exphcit dependence 
has been suppressed for brevity. We use the definitions of 6 and K as they are given 
in ref. these definitions differ from the ones used in refs. |^, ^, ^ by a sign in B and 
a factor of 2tt in K. 
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Jo 


Jn (n^O) 


1> A> B 
A>1> B 
A> B>1 

1> B> A 
B >1> A 
B> A>1 


2mC 

2m{\nA + C) 
2m{\uA-\nB + C) 

2mC 

2m{-\YiB + C) 
27ri(ln^-ln5 + C) 


27ri(E" - 1) 


27ri(E" - A")/n 
2m\l - A") 




Table 11. The integrals J, 



26 



Figure captions 



Fig. 1. The six vertices and their energy assignments. 



Fig. 2. Phase diagram of the six-vertex model in the {h, v) plane for the case A > 1. 
(Specific parameter values are S — 1.5, e = 0.3, and /c^T = 1.1.) The bold lines 
show the phase boundaries. The dashed lines show the division of the incompletely 
polarized part of the phase diagram into the six regions which are discussed in the 
text. 



Fig. 3. Free energy surface F{h, v) and phase diagram for the six-vertex model in the 
{h, v) plane for the case A > 1. (Specific parameter values are S = 1.1, e = 0.45, and 
ksT = 0.35.) The bold lines on both surface and base plane show the phase bound- 
aries. For the completely polarized phases, the values of the polarizations {x, y) are 
noted. All boundaries are second-order (specifically, Pokrovsky-Talapov), except 
for the one between the {x,y) = (1,1) and (—1,-1) phases, which is first-order. 
The filled squares at each end of the first-order line are the conical points at 
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Fig. 4. A section of the phase diagram shown in the base plane of Fig. |^. The sohd 
hnes again denote phase boundaries and the filled square denotes the conical point. 
The dashed and dotted lines are some contours of constant a and 6, respectively. 
From top to bottom, the dashed lines show a ^ —0.07, —0.20, —0.51, and —1.02. 
From left to right, the dotted lines show h = {y + 0o)/2 ~ 2.59 (which gives y = —x 
and h = —v), b ^ 2.72, 6 = 0o ~ 2.87 (which gives x = 0), and b ^ 3.13. Note that 
the phase boundaries themselves correspond to the limit a — 0, with the value of b 
determining the location along these boundaries. The conical point corresponds to 
either a —it or 6 — » oo. In the former limit, the angle of approach to the conical 
point is determined by the value of b, while in the latter limit, the approach is tangent 
to the phase boundary for all values of a. 



Fig. 5. Free energy surface F{x,y) and phase diagram for the six-vertex model in 
the {x, y) plane, for the same parameter values as in Fig. ^. The two coexistence 
regions correspond to the two conical points {h,v) = ±(— ^) at each end of 
the first-order line in Fig. The regions are bounded by the lines l± and y = x 
(shown as bold lines). The line y = x is also a coexistence region, corresponding to 
the entire first-order line in Fig. ^. [Although it is not obvious from this perspective, 
F{x,y) is constant along y = x.] 



Fig. 6. Phase diagram of the six-vertex model in the {h, v) plane for the case A = 1. 
(Specific parameter values are 6 = 1.1, e = 0.45, and ksT = ksTc ~ 0.5913.) The 
bold lines show the phase boundaries. For the completely polarized phases, the values 
of the polarizations (x, y) are noted. The phase boundaries are all second-order, since 
the first-order boundary between the {x,y) = (1, 1) and (—1, —1) phases has shrunk 
to zero length. The two conical points have thus coalesced into a single point at 
ih,v) = iO,0). 



Fig. 7. Contour of integration in the complex 2;-plane for the integrals J„ of Eq. (p.l| ). 
A branch cut runs along the real axis between A and B. The labels and 27r above 
and below this branch cut give the value of the imaginary part of the logarithm. For 
the case n = 0, there is also a pole at z = 0. 

Fig. 8. Contour of integration in the complex 2;-plane for the integrals Jn of Eq. ( |C.1| ) 
for the case B < A < 1. A branch cut runs along the real axis between A and B. 
The labels ivr above and below this branch cut give the value of the imaginary part 
of the logarithm. For the case n = 0, there is also a pole at ^ = 0. 
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